Solving frustration-free spin systems 
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We identify a large class of quantum many-body systems that can be solved exactly: natural frustration-free 
spin-1/2 nearest-neighbor Hamiltonians on arbitrary lattices. We show that the entire ground state manifold of 
such models can be found exactly by a tensor network of isometries acting on a space locally isomorphic to the 
symmetric subspace. Thus, for this wide class of models real-space renormalization can be made exact. Our 
findings also imply that every such frustration-free spin model satisfies an area law for the entanglement entropy 
of the ground state, establishing a novel large class of models for which an area law is known. Finally, we show 
that our approach gives rise to an ansatz class useful for the simulation of almost frustration-free models in a 
simple fashion, outperforming mean field theory. 



Understanding the physics of quantum many-body systems 
is a central goal of modern physics, as they can exhibit ex- 
otic phenomena with no parallel in classical physics, includ- 
ing topological effects and quantum phase transitions at zero 
temperature. However, the very source of their rich physics 
also leads to a major roadblock in their study: the Hilbert 
space dimension of these systems scale exponentially with the 
number of particles. This means that brute-force numerical 
techniques fail even for systems of only a handful of particles. 

A key insight in the study of local quantum many -body sys- 
tems is that naturally occurring states only occupy a small 
subspace of the Hilbert space which in principle is available 
to them. Specifically, it has been realized that ground, ther- 
mal, and dynamically evolving states are only weakly en- 
tangled: the entanglement entropy satisfies what is referred 
to as an "area law" HHH]. This insight is the basis of the 
density-matrix renormalization group approach, and higher- 
dimensional analogues [4]. So successful are these methods in 
practice that one is tempted to boldly conjecture that all physi- 
cally relevant systems will soon be tractable to one or another 
of the numerical tools we have to hand. Recent important 
results give cause for caution, showing that general numeri- 
cal methods cannot well-approximate the physics of an arbi- 
trary local quantum system, even in ID: these include local 
glassy models where approximating the ground state energy 
is QMA-complete |5|, suggesting that such systems would be 
intractable even for quantum computers. However, these re- 
sults do not give much reason for practical concern thus far, 
as at least in ID they rely on rather baroque constructions (in- 
volving very large local dimensions). 

In this work, we approach the issue of the complexity or 
"hardness" of finding ground states from the other direction: 
We establish a large class of models for which the task of find- 
ing the ground state is easy, in that the ground state mani- 
fold can be described exactly and efficiently. This is the class 
of all natural frustration free spin-l/2 models with nearest- 
neighbor interaction on general lattices. (By the qualifier 
"natural", we mean that all two-spin interaction terms have 
excited states which are entangled, which might be taken as 
implicit in the semantics of an "interaction term".) Adopt- 
ing and extending ideas of Ref. [6] on QUANTUM 2-SAT and 



going beyond ID models as in Ref. [7|, we find that the com- 
plete ground state manifold of such Hamiltonians can be con- 
structed by considering the ground spaces of each interaction 
term in turn, and obtaining a reduction to the symmetric sub- 
space of a smaller system. In doing so, we find that the re- 
sulting ground state manifold can be efficiently grasped in 
terms of tree-tensor networks. We discuss how this allows 
expectation values of local observables to be computed effi- 
ciently. What is more, the ground states satisfy an area law. 
Physically, we can view this work as describing a large class 
of models for which an instance of real-space renormaliza- 
tion provides an exact solution to the true genuine quantum 
many-body model. Finally, we see how this construction — a 
tree tensor network with a symmetric subspace as an input — 
can serve as an ansatz class to simulate almost frustration-free 
models which are in a sense "close" to frustration-free mod- 
els. In this way, one can outperform mean-field approaches in 
a very simple fashion. 

In our analysis, we allow for nearest-neighbor Hamiltoni- 
ans on arbitrary lattices. This could be a cubic lattice of some 
dimension, or more generally any graph, the vertex set of 
which we denote by V. On this lattice, the spin Hamiltonian 
H is represented as 
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for terms h a ,b acting on pairs of spins {a, b} C V. By rescal- 
ing, we may without loss of generality require that the ground 
energy of each term h a j, is zero. 

The ground state manifold M of such a Hamiltonian may 
be degenerate: we identify the ground state p with the maxi- 
mal mixture over M . (This is a pure state only if H is non- 
degenerate.) We describe properties of the ground state p, 
and more generally the manifold M, given H as in Eq. ([TJ. 
The Hamiltonian H is frustration-free (or unfrustrated) if the 
ground state vectors |$) £ M are ground states of individual 
coupling terms, that is if h a j,\^) = holds for all h a ^ and all 
|<&) £ M; we say otherwise that H is frustrated. 

We will call a spin Hamiltonian H natural if it contains no 
isolated subsystems, and each interaction term h a ,b (consid- 
ered as an operator on C 2 ® C 2 ) has at least one entangled 



excited state (i.e., an entangled state orthogonal to the ground 
state manifold of h a ^). In what follows we will consider only 
such "natural" Hamiltonians. 

Frustration-free spin Hamiltonians. Our main results con- 
cern the class of (natural) frustration-free spin Hamiltoni- 
ans H. We show that the ground state manifold M of a 
frustration-free spin Hamiltonian on N spins has dimension 
at most N + 1. What is more, it is the image of a space of low 
Schmidt measure [ 8 1 under a tree-tensor network. We proceed 
by extending the work of Ref. [6| for QUANTUM 2- SAT. 

First, we describe the needed components of Ref. [6|, in the 
language of frustration-free models. Consider a Hamiltonian 
Hij containing terms h UjV of rank 2 or 3. If Hu is frustration- 
free, the reduced state p U:V of any state vector |$) G ker(_ff) 
is in the kernel of h u>v ; we may then consider a subspace 
S u , v C H u ® H v of dimension 2 which contains supp(p Ui „). 
By defining an isometry 
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we can reduce to a Hamiltonian on fewer spins: we let 



{a,b} 

Such a spin Hamiltonian H' v is a sum of two-spin interac- 
tions (and possibly single-spin terms) of the form h' ab = 

Ru-.uvha.bRu-.uv (If h u ,v has rank 3, then h' uv is a non-zero 
single spin operator acting on u alone.) If H contains non- 
zero terms h au and h av , we obtain two non-zero contribu- 
tions fl au — ti u:uv ila,u-L*'u:uv and fl av — ^ u:uv tla,v^u:uv In 

the Hamiltonian H\ each of which act on {u, a}. The sum 
gives a combined term h' au — h' au + h' av in H' v , possibly of 
higher rank than either h' a u or h' a v 11151 . If the new Hamil- 
tonian H'jj contains terms of rank 2 or 3, we may perform 
another such reduction, and so on. This reduction procedure 
has the following features: 

a. Preservation of the kernel dimension. By construc- 
tion, we have ker(Hu) = R u:uv ker(iJ^). Thus, the kernels 
of Hjj and H[j have the same dimension. In particular, if H' [f 
has any terms of full rank acting either on one or two spins, 
then dimkei (Hu) = 0, in which case Hjj is frustrated. If 
no full-rank terms are produced, each reduction leads to an 
operator acting on fewer spins, until we obtain a Hamiltonian 
having only terms of rank 1, 

b. Arbitrariness of reduction order. Because the dimen- 
sion of the kernel is preserved by these reductions, we may 
continue to perform such reductions until we obtain a Hamil- 
tonian which either (/) contains only terms of rank 1, or 
(//) contains a full-rank term. The latter cannot occur in the 
reduction of a frustration-free Hamiltonian; and we discuss 
below the analysis for Hamiltonians having only rank-1 terms. 
Thus, we may choose any convenient reduction sequence. 

The above features allow us to reduce to the special case 
of a Hamiltonian H v (acting on a system V') which has only 
interaction terms of rank 1 . Each two-spin Hamiltonian term 
h' a b ~ \Pa,b)(Pa,b\ ma Y be regarded as imposing constraints 
on the corresponding two-spin marginals p a j, of states |$) € 



ker(H): we aim to obtain additional constraints on pairs of 
spins u, v € V by combining the known constraints. To this 
end, Ref. ]6) shows that a state |$) which is in the kernel of 
two functionals (/3 a ,b\,{Pb,c\ i s a l so m the kernel of 



(P'a,c\ = ({P;b\®{Pb,e\)( 1 ®\*~)® 1 ) 
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acting on the spins {a,c}, where | X I'~) is the two-spin an- 
tisymmetric state vector. For each such "induced" con- 
straint (j3' u J on spins {u,v}, we may add the term h UtV = 
\P'u vriP'u v\ t° Hjj, resulting in a Hamiltonian Hu which has 
the same kernel as H[j. (If H'y contains a term h' uv gk h u ,v> 
it can be subsumed into a term h u>v with rank at least 2, in 
which case we apply a reduction R u:uv as above.) One may 
induce further constraints from the terms of H, until we arrive 
at a "complete homogeneous" Hamiltonian H c , having only 
terms of rank 1, for which the constraints (/3 U „ | are closed (up 
to scalars) under the constraint-induction procedure of Eq. (01. 
Ref. [6j shows that such a Hamiltonian H c , acting on at 
least one spin and lacking single-spin operators [17|, has a 
ground space containing product states. Thus, the above re- 
marks essentially recap the following result: 

Observation 1 (Ref. |6|): There is an efficient algorithm to 
determine whether a spin Hamiltonian is frustration-free. 

We now extend the above results, to obtain a strong char- 
acterization of ground state manifolds for natural frustration- 
free systems. We note the following three additional features 
of the isometric contraction scheme above: 

c. Tree-tensor construction. The complete network of 
isometric reductions T represents a tree-tensor network, a spe- 
cial case of the MERA ansatz 1 1 1 1 which is related to real- 
space renormalization. Each isometry R u -. uv has one free in- 
put tensor index and two free output indices, and the sequen- 
tial nature of the reduction ensures that the network is a di- 
rected and acyclic. Thus, any spin v introduced by an isom- 
etry Ru-uv is a "daughter spin" of a unique parent u, leading 
to a tree-like structure on the tensor network T. Note how- 
ever that T has free input indices, corresponding to the roots 
of each tree: by construction, the ground space of Hu is the 
image T|*) of states |*) e ker(H c ). 

d. Preservation of natural Hamiltonians. Importantly, 
the isometric reductions above preserve the class of natu- 
ral frustration-free spin Hamiltonians: that is, the mapping 
h a ,u >-> Rt:uvha,uRu:uv does not decrease the rank of the in- 
teraction on {a, u], and does not map the orthocomplement 
of the kernel to a space of the form | <fi) ® C 2 for any | <j>). 

e. Reduction to the symmetric subspace. As a conse- 
quence of the previous feature, we may use the isometric re- 
ductions to map any "natural" frustration-free Hamiltonian to 
a "complete homogeneous" Hamiltonian H c on a system V c , 
in which the non-zero terms h a ,b = \Pa,b)(/3 a ,b\ are supported 
on entangled states |/3 a ,b)- We show that the kernel of such 
a Hamiltonian has small dimension, and is spanned entirely 
by product states. For an arbitrary spin a 6 V c , we may let 
L a = 1 and define a family of operators L v satisfying 
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for spins v E V c and operators (j3 a , v \. One then finds that 
C = <S>v£v L v is a linear isomorphism (not necessarily 
an isometry) from the subspace Symm('Hf 1 c ) of symmetric 
states to the ground space of H c . (This isomorphism is also 
noted in Ref. [12 Section III A] for generic Hamiltonians with 
rank-1 interactions.) As Syram^Hf ) i s spanned by uniform 
superpositions \Wk) of standard basis states having Hamming 
weight O < n c = \V C \, we have dimSymm(^f nc ) 
n c + 1. This subspace may also be spanned by product state 



Observation 3: The ground states of any frustration-free spin 
Hamiltonian Hjj on a lattice obey an entanglement area law. 

For any contiguous subsystem A containing a spins, we 
may reduce the Hamiltonian H\j acting internally on A — 
by a tree-tensor isometry Ta acting on A alone, and inducing 
constraints as described by Eq. (RJ — to obtain a complete 

(A) 



_i homogeneous Hamiltonian He , acting on at most a spins. 

(A) 

The ground space of He has dimension at most a + 1; as 

. , \a nc )® n " for any set of n c + 1 pairwise the ground space of Hjj ' is an isometric image of that of 

independent state vectors \ aj ) e H 2 - Thus, any complete #) the same is true of H { T f ] . As H„ is frustration-free, 
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homogeneous (natural) Hamiltonian H c has a ground space 
spanned by a family of vectors 



I*,- 



L, : 
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v£V c 



for some |ao), ■ ■ ■ , \<Xn c ) € H.2 as described. 

Coupled with the tree-tensor structure of the isometric re- 
ductions, this characterization of the ground space of (natural) 
complete homogeneous spin Hamiltonians has the following 
consequences for (natural) frustration-free Hamiltonians: 

Observation 2: For any frustration-free spin Hamiltonian 
Hy, any constant k, and for k-local operators A, (A) can 
be efficiently computed with respect to ground states of H. 

Let H c be a complete homogeneous Hamiltonian obtained 
by isometric reduction of an unfrustrated Hamiltonian Hjj, 
and let n c be the number of spins on which H c acts. Con- 
sider a fc-local operator A. As ker(_ff c ) is spanned by some 
collection of product vectors |* ) = C|a )®™ c , . . . , |^„J = 
C|a„ c )®™ c as in Eq. d6j, we can efficiently compute the re- 
striction of A to ker(i/ c ) by evaluating the matrix 



W(A) = J2 liX*;l4**> 

j,k=0 



(7) 



followed by a suitable transformation. Specifically, consider 
the operator B = W(l); we have B = UAW for some U 
unitary and A positive and diagonal. We find 

A- 1/a J7t£uX*i| = £lj'X*il. (8) 

for some orthonormal basis |$o), . . . , \®n c ) of ker(H c ); thus, 
the restriction of A to ker(iJ c ) with respect to the basis of 
states |$j) may be computed as 



A =A- 1/2 U' i W(A)UA- 1 



/2 



(9) 



(For A consisting of a single fc-spin term, the inner products 
of Eq. |7]l consist of a product of constant-dimensional inner 
products; for A a sum of multiple terms we extend linearly.) 
Let T : ker(H c ) —> kei(Hjj) be the network of isometric re- 
ductions: then, by considering operators A = T^ AT, we may 
compute the restriction A of such operators A to the ground 
space of Hy. We may then efficiently compute estimates us- 
ing such (polynomial-size) matrices. 



(A) 

any ground state of Hjj is also a ground state of H\j ; it then 
follows that the Schmidt measure [8] of \ty) with respect to 
the bipartition V — A U (V \ A) is at most log (a + 1). For a 
lattice with a well-defined dimension, we may obtain an area 
law for arbitrary subsystems of the spin lattice by summing 
this logarithmic bound over the number of distinct connected 
components. 

"Almost" frustration-free Hamiltonians. We now leave the 
rigorous exact setting and turn to the physically interesting ob- 
servation that the above ground state manifold can be used to 
grasp approximately frustration-free models. Indeed, Obser- 
vation [2] in particular suggests a variational approach to esti- 
mating ground energies for Hamiltonians H for which 



H = H,, 



■XHj, 



(10) 



for some A <C 1, where H\j is frustration-free but the Hamil- 
tonian H itself is frustrated. For such Hamiltonians H, if 
no phase transition is encountered, the eigenvalues and eigen- 
states may differ little from those of the unfrustrated Hamil- 
tonian Hjj (for related bounds on eigenvalues, see, e.g., Ref. 
lfT3l ). If the lowest k eigenvectors (for some suitable 1 < 
k < n + 1) have a sufficiently high overlap with the lowest k 
eigenvectors of Hjj, we may approximately sample from the 
low-energy eigenvectors of H by restricting to the kernel of 
Hjj, using the efficient algorithm above. In particular, as this 
procedure is variational, estimates obtained in this way for the 
ground state energy of H are guaranteed to be upper bounds. 

Observation 4: The ground state manifolds of frustration-free 
Hamiltonians serve as an ansatz class for almost frustration- 
free models. 

We may also consider additional improvements to this es- 
timation ansatz, in which less information about the "frus- 
trating" component XHp of the Hamiltonian is lost than in 
the tree-tensor renormalization procedure for the frustration- 
free model Hjj. To obtain better estimates — and to extend 
these techniques to the case where A may be significantly large 
— we may consider (a) partial reductions by tree-tensor net- 
works, having many free input spins; (b) isometries depend- 
ing on the terms of the perturbed Hamiltonian H (rather than 
those of Hjj); and (c) using additional variational approaches. 
We describe these ideas below. 

In the case that the ground state manifold of Hjj is con- 
tained in a subspace K, which is spanned by product states 
and has "small" dimension (i.e., polynomial rather than ex- 
ponential in the system size, as in the case of the symmetric 
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FIG. 1. Left: Ground-state energy for XXZ-model on a trigonal 
lattice on a 3 x 3 torus, hij = — XiXj — Y t Yj — (1 — X)ZiZj, 
by symmetric subspace estimate compared to product state ansatz 
and exact diagonalization. The inset shows the same model on a 
6x6 torus where exact solution is not feasible and, therefore, is 
replaced by an Anderson lower bound. Right: Magnetization in 
z-direction for Ising model in a transverse field on a 4 x 4 torus, 
hij — —ZiZj, hi = —XXi, by symmetric subspace estimate com- 
pared to Gutzwiller mean-field approximation (product state ansatz) 
and exact diagonalization. 



subspace for ferromagnetic Ising or XXX models), we may 
indeed forgo isometric reductions entirely, and estimate the 
ground energy of H by considering the restriction of H to 
K, using the techniques described for Observation [2j In con- 
trast to a full tree-tensor contraction, this has the advantage of 
yielding exact results for the frustration-free case A = while 
retaining more information about the "frustrating" component 
H p . The resulting estimate for the ground state energy will be 
a linear function of A, whose value and first derivative with re- 
spect to A agree with that of the exact ground energy at A = 0. 
For "small" A and modest system sizes, this may yield a good 
estimate of the ground state energy of H: see Fig.[T]|14]. 

To obtain estimates which account for spatially decaying 
correlations, we may perform a partial tree-tensor reduction 
with a small number of contraction layers, and sample with 
respect to a subspace K, as described above. In each layer, we 
may fix a collection of (non-intersecting) adjacent site pairs 
{o, b} to contract, and for each such pair {a, b} apply some 
isometric contraction as described in Eq. ([3]). However, rather 
than apply the reductions which would be suggested by the 
frustration-free Hamiltonian Hjj, we may use isometries 



}a:ab = |< 6 )<0| + |< b }(l| 



(ID 



where \ip® b ), \ifj^ b ) are the lowest energy eigenvectors of r\ a &, 
where r\ a j is the partial trace of H with respect to all sites 



other than a and b. Given a tensor network T\ consisting of 
a product of such two-site operators Q a -.ab, we may then con- 
sider the spin model given by H' = T{HTi, and estimate 
the ground energy of H' with respect to a low-dimension sub- 
space or another isometric contraction; this yields an upper 
bound on the ground energy of the original Hamiltonian H. 

As expected, the symmetric subspace estimate works best 
for almost-unfrustrated systems. The region where the results 
are accurate is larger in small and medium size systems, where 
the improvement compared to mean-field is most significant: 
see Fig.fTI For the XXZ-model, applying isometric reductions 
on disjoint pairs of nearest neighbors improves the approxi- 
mation dramatically (versus symmetric states) for "strongly" 
frustrated interactions, as they allow for the selection of the 
most important two-site subspaces fl31 . 

The approximation technique presented here can also be 
combined with any variational method. One simple idea is 
to apply single-site unitaries to the Hamiltonian before per- 
forming the restriction to the symmetric subspace. One may 
then vary the (different) unitaries applied on the sites to mini- 
mize the energy, in order to find a better approximation to the 
ground state. This can lead to a notable improvement, espe- 
cially when antiferromagnetic effects are non-negligible. By 
construction, this approach performs at least as well as a prod- 
uct state ansatz to which it reduces for A->oo [17|. 

Summary. In this work, we have introduced a class of spin 
models that can be completely solved: The entire ground state 
manifold can be explicitly given and parametrized by an en- 
tire symmetric subspace under a tensor network, so a "tensor 
network with an input". This class of models is hence ex- 
pected to constitute a rich playground of exploring ideas on 
quantum lattice models, complementing work that exempli- 
fies how computationally difficult it can be to even approxi- 
mate ground state energies. In a sense, the considered models 
can also be viewed as the parent Hamiltonian of the network, 
in a converse approach taken for tree tensor networks in Ref. 
lfT8ll . We always find area laws for the entanglement entropy, 
hence establishing a new large class of models beyond free 
systems for which an area law can be proven to hold. It is the 
hope that this work stimulates further research on models for 
which tensor networks not only arise as computational but as 
essentially analytical tools. 
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